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1 Introduction 

There is no doubt that the Standard Model (SM) 
describes physical processes at the highest avail- 
able energies with very good precision. Yet it con- 
tains some 18 independent parameters which re- 
main completely undetermined and await deriva- 
tion from a more fundamental theory. 13 of these 
parameters are related to the fermion mass sec- 
tor which is clearly distinguished by an amazingly 
simple regularity in the hierarchy of masses and 
mixing angles of the three fermionic families. 

This observed pattern may provide a clue to a 
more fundamental theory at some higher scale M 
(M ~ Mpi anck or MstHng or Mqut)- Our main 
hypothesis is that only a small set of effective mass 
operators below M dominates in the Yukawa ma- 
trices, subsequently leading to the fermion mass 
matrices at the scale of electroweak symmetry 
breaking. These effective operators are of utmost 
interest. They can lead to the reconstruction of 
the full effective field theory at the scale M and 
provide the matching conditions for a fundamental 
string theory, believed to be the ultimate quan- 
tum theory incorporating gravity together with 
the gauge interactions of quarks and leptons. 

SO(10) supersymmetric theories are excellent 
candidates for such a theory below the string scale. 
They maintain the successful prediction for gauge 
coupling unification and provide a powerful frame- 
work for predictive theories of fermion masses and 
mixing angles. This is because all the fermions of 
a single family are contained in the 16 dimensional 
representation of SO(10) - thus fermion mass ma- 
trices are related by symmetry. In the most pre- 
dictive theories, the ratio of Higgs vevs - tan/3 - 
is large, and the top quark is naturally heavy as 
found experimentally. However, as a consequence 



of large tan/3 there are potentially large super- 
symmetric weak scale threshold corrections which 
could play an important role in fitting the fermion 
masses and mixings. Thus a self-consistent analy- 
sis necessarily includes the dimensionful soft SUSY 
breaking parameters, in addition to the dimension- 
less gauge and Yukawa couplings. 

2 Connecting GUT Scale Physics with 
Low Energy Observables 

Here we present the results of such a complete 
top-bottom analysis. It starts at the GUT scale 
Mq, which is a free parameter itself, with uni- 
fied gauge coupling ac, n y free parameters enter- 
ing the Yukawa matrices and with five universal 
soft SUSY breaking parameters fx, mo, M1/2, Aq 
and B[i. In addition, we introduce €3 as a one 
loop GUT threshold correction to a s (Mo), and 
non-universal Higgs masses mn d and 171h„ ■ The 
dimensionless (dimensionful) couplings are run 
down to the Z scale using two (one) loop renor- 
malization group equations (RGEs) of MSSM. Wc 
have checked at a few selected points that the cor- 
rections to our results obtained by using two loop 
RGEs for dimensionful couplings are, in fact, in- 
significant. At the Z scale we match the MSSM di- 
rectly L!l to the SU(3) C x U(1)em, thus leaving out 
the SM as an effective theory on our way down 
to the experimentally measured low energy data 
listed in table [l]. At the Z scale we use the tree 
level conditions for electroweak symmetry break- 
ing to fix v and tan/3. Then, still within the 
MSSM, we calculate the amplitude for the pro- 
cess b — > sj , one loop corrected W and Z masses, 

"Clearly, n y is model dependent. 

b By definition, Mq is to be understood as the scale 
where ct\ and 012 are the same, and equal to ag . 
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Table 1 : Experimental values of the twenty low energy ob- 
scrvablcs entering the \ 2 function, together with the cor- 
responding standard deviations. 



Observable 


Central value 


a 


M z 


91.186 


0.46 


M w 


80.356 


0.40 




0.11664 x 10~ 4 


0.0012 x 10~ 4 


a EM 


137.04 


0.69 


a s (M z ) 


0.118 


0.005 


Pnew 


-0.6 x 10" 3 


2.6 x 10~ 3 


M t 


175. 


6. 


m b (M b ) 


4.26 


0.11 


M b - M c 


3.4 


0.2 


m s 


180. 


50. 


m d /m s 

Q- 2 


0.05 


0.015 


0.00203 


0.00020 


M T 


1.777 


0.0089 


M„ 


105.66 


0.53 


M e 


0.5110 


0.0026 


V us 


0.2205 


0.0026 


Vcb 


0.0392 


0.003 


V u b/Vcb 


0.08 


0.02 


B K 


0.8 


0.1 


BR(b -> s 7 ) 


0.232 x 10~ 3 


0.092 x 10~ 3 



corrections to the p parameter from new physics 
outside the SM, and one loop corrected G M , where 
we leave out the SUSY vertex and box contribu- 
tions to the Ar parameter. When crossing the Z 
scale, we compute the complete one loop threshold 
corrections to a s and aEM, whereas only those one 
loop threshold corrections to the fermion masses 
and mixings enhanced by tan/3 are computed. □ 
Masses and couplings which require further run- 
ning are renormalized down to their apropriate 
scales using three loop QCD and one loop QED 
RGEs. The branching ratio BR(b — ► sj ) is renor- 
malized down to the scale Mb using the leading 
log approximation. Note that some observables 
entering the x 2 function (see table 0) are known 
so well that we have to assign a theoretical, in- 
stead of experimental, error as their standard de- 
viation. This is true for Mz, Mw, olem, G> and 
the charged lepton masses. We estimated conser- 
vatively our theoretical error to be 0.5% based on 
the uncertainties from higher order perturbation 
theory and from the performance of our numeri- 
cal analysis. The error on G p is estimated to be 
within 1% due to the fact that in addition to the 
uncertainties mentioned above we neglect SUSY 



vertex and box corrections to Ar. Also note that 
€k, the observable of CP violation, has been re- 
placed by a less precisely known hadronic matrix 
element Bx- Thus our theoretical value of -Br- is 
defined as that value needed to agree with for 
a set of fermion masses and mixing angles de- 
rived from the GUT scale. Finally note, that the 
light quark masses are replaced by their ratios (at 
the scale lGeV), and m c (M c ) by the difference 
Mb — M c , since the latter quantities are known 
to better accuracy. Q~ 2 — (m 2 , — m\)/m 2 s is the 
Kaplan-Manohar-Leutwyler ellipse parameter. 

In addition, the % 2 function is increased signif- 
icantly by a special penalty whenever a sparticle 
mass goes below its current experimental limit. 

3 Model 4 of ADHRS 

We have analyzed several models of fermion 
masses. First, we. have studied the best working 
model of ADHRSB, model 4. The model is defined 
by four effective operators at the GUT scale (A's 
stand for adjoint states and S for a singlet state): 
33 = 16 3 1O16 3 , 023 = 16 2 10 4* 16 3 , 

012 = 16i (|) 3 10 (|) 3 16 2 and by one of 
six possible 022 operators (they all give the same 
0:1:3 Clebsch relation between up quarks, down 
quarks and charged leptons responsible for the 
Georgi-Jarlskog relation). Each operator enters 
with its own complex coefficient, and we can rotate 
away three independent phases. Hence in model 4 
n. y =5; the number of initial parameters for the \ 2 
analysis is 15 and we are left with 5 degrees of free- 
dom. Fig. [l] shows the contour lines of constant x 2 
in the mo — Mu2 plane for different values of /i. 
No substantial improvement of the performance 
of model 4 can be achieved by neglecting one out 
of the twenty low energy data given in table [y. □ 
On the other hand we found that a significant im- 
provement is possible by adding one new opera- 
tor, contributing to the 13 and 31 elements of the 
Yukawa matrices. 

4 Model 4c 

Next, we analyzed two models derived from the 
complete SO(10) SUSY GUTs discussed recently 
by Lucas and Rabp. The models were constructed 
as simple extensions of model 4. Different label 
(a,b,...f) refers to the different possible 22 oper- 
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Figure la: M=80GeV 
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Figure lc: M=240GeV 
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Figure 1: Model 4 global analysis results for different values 
of fi(Mz)- Solid (dashed, double-dash-dotted) lines repre- 
sent contour lines of constant \ 2 = 15 ( 14, 13 ) / 5dof. 
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Figure 2: Model 4c global analysis results for different val- 
ues of ii(Mz)- Solid (double-dash-dotted, dotted) lines rep- 
resent contour lines of constant x 2 = 6 ( 3, 1 ) / 3dof. 



ators.j] In the extension to a complete GUT the 
different 22 operators lead to inequivalent theories 
due to different U(l) charge assignments. When 
one demands "naturalness" , i.e. includes all terms 
in the superpotential consistent with the symme- 
tries of the theory one finds only one new oper- 
ator (0i3) for models 4a and 4c. Model 4b, on 
the other hand, has no new operator and thus is 
equivalent to model 4 (above). The 22 and 13 op- 
erators of model 4c are O22 = I62 ^ 10 4^- I62 and 
Oi3 = I61 (|) 3 10 4f 16 3 . With the 13 operator 
n y =7, which implies 3 degrees of freedom. The 
results of the global analysis are given in Fig. ||, 
and in table || for a selected point II marked on 
Fig. ||a. Model 4a is defined by different 22 and 



c Note, models d, e and f have the second family I62 
coupled directly to 10i and a heavy 16. If this coupling 
is as large as the third generation Yukawa coupling, then 
we would obtain excessively large flavor changing neutral 
current processes, such as /i — » cy. Thus these models were 
not considered in ref. [5]. 



13 operators and gives the fits with the best x 2 — 4- 
6 in most of the SUSY parameter space. (It yields 
also x 2 — 3 j but only for the corner in the SUSY 
parameter space with large mo, M\/2 and fi.) □ 

5 Discussion and Conclusions 

The results of our analysis, as well as the whole 
project presented here, may be understood from 
two different perspectives. 

The emphasis of this work has been put on 
the study of the origin of fermion mass matri- 
ces, within the context of minimal SO(10) SUSY 
GUTs. The best working model of ADHRS, with 
just four effective operators in the Yukawa sec- 
tor, can be assigned a confidence level of only 3- 
4%, in spite of the fact that each observable is 
within 2(7 of its experimental value. The addition 
of a 13 operator may improve the performance of 
the model. Substantial improvement however, is 
not automatic, as was mentioned in the example 
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Table 2: Model 4c - Results at Point II. Initial parameters: 
1/og = 24.36, Mo = 3.1710 16 GeV, e 3 = -4.89% , A = 
0.807, B = 5.44-10" 2 , C = 1.15-10" 4 , D = 4.94-10" 4 , <5 = 
5.71, E = 1.3M0- 2 , 3> = 1.04, fi = 80GeV, m = 700GeV, 
M 1/2 = 240GeV, m Hd /m = 1.42, m Hu /m = 1.24, A = 
458.35GeV, Bp. = 120.66GeV 2 . The last column displays 
SUSY threshold corrections in %. 



Observable 


Computed 


Partial \ 2 


S.t.c. 


M z 


91.12 


0.02 




M w 


80.38 


<0.02 




G> 


1.166-10" 5 


<0.02 




a EM 


137.0 


<0.02 


1.43 


a s (M z ) 


0.1151 


0.34 


12.78 


Pnew 


+ 1.87-10- 4 


0.09 




M t 


175.7 


<0.02 


0.74 


m b (M b ) 


4.287 


0.06 


5.43 


M b - M c 


3.440 


0.04 


7.56 


m s 


189.0 


0.03 


3.68 


m d /m s 

Q- 2 


0.0502 


<0.02 


0.00 


0.00204 


<0.02 


1.78 


M T 


1.776 


<0.02 


-2.08 


M„ 


105.7 


<0.02 


-1.50 


M e 


0.5110 


<0.02 


-1.50 


v us 


0.2205 


<0.02 


0.00 


v cb 


0.0400 


0.07 


1.58 


v ub /v cb 


0.0772 


<0.02 


0.00 


B K 


0.8140 


<0.02 


-3.18 


BR(& -> s 7 ) 


2.382-10- 4 


<0.02 




TOTAL x 2 




0.7306 





of model 4a. Nevertheless, we showed that model 
4c provides an excellent fit to all 20 observables 
of table [lj with the confidence level better than 
68% in a large region of the experimentally al- 
lowed SUSY parameter space. Note that the best 
fits extend to the region with very large mo where 
the effect of the SUSY threshold corrections to 
fcrmion masses and mixings is suppressed by large 
squark and slepton masses. As a result, in this re- 
gion the effective number of degrees of freedom in 
the fermionic sector is actually larger than 3, since 
there are 7 parameters (A, B, C, D, 5, E, $) in the 
Yukawa matrices determining the 13 low energy 
masses and mixings of the fermions. This means 
that the Yukawa sector of model 4c does actually 
a much better job than appears at first glance. 
Whether or not this particular model is close to 
the path Nature has chosen remains to be seen. 
One important test will be via the CP violating 
decays of the B meson. Model 4c predicts a value 



for sin2a ~ 0.95 which is insensitive to the SUSY 
breaking parameters, whereas in the SM the value 
of sin2a is unrestricted]. Another important test 
may come from nucleon decay rates 0. 

When one disregards the origin of the Yukawa 
matrices at the GUT scale, this work can also be 
viewed as an MSSM global fit in the large tan/3 
regime. From this perspective, a special feature 
of our approach is that we are getting complete 
3x3 fermion and sfermion mass matrices at the 
SUSY breaking scale, instead of just the leading 
33 elements. Viewed as MSSM global fits, our 
results suggest that there is no narrow, strongly 
preferred region in the SUSY parameter space. 
Hence, at present one cannot make strong con- 
clusions (with tan/3 large) about the masses of 
sfermions, which leaves open manji channels for 
Tevatron and LEPII experiments. □ It is also in- 
teresting to note that the best fits at a number 
of points give a low pseudoscalar Higgs mass and 
low valuesa of a s (Mz) leaving the door open for 
a natural explanation of the 1-2ct increase in the 
partial width Z — > bb, without any tuning of the 
initial parameters to obtain this effect. Finally, 
by keeping the complete 3x3 mass matrices for 
both fermions and sfermions one can study flavor 
dependent processes in a theory which fits the low 
energy data; for example, rare B and K decays, 
B — B mixing, lepton flavor violating processes or 
even Z — > bb. 

At the present time, when direct evidence for 
physics beyond the SM evades experimental ob- 
servations, global analysis serves as the best test 
of new physics. Moreover, our results show that 
minimal SUSY SO(10) models remain among the 
best candidates for the theory of nature beyond 
the Standard Model. 
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